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the variability properties of a large homogeneous sample of SDSS quasars. We search for an excess 
of characteristic variations of quasar brightness that are associated with string lensing and use the 
observed distribution of this variation to constrain the density of strings. The limits obtained do not 
invoke any clustering of strings, apply to both open segments and closed loops of strings, usefully 
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and below. Further progress in this direction will depend on better understanding of quasar intrinsic 
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I. INTRODUCTION 

Linear topological defects arise naturally during phase 
transitions in diverse areas of physics. Various processes in 
the early universe could also produce such defects, which are 
called cosmic strings It is often assumed that phase 

transitions lead to formation of strings with a characteris- 
tic tension of order the squared energy scale of the string- 
producing theory (in Planck units). After formation, the 
strings build up an intricate network, combined from open 
segments of the horizon scale and a multitude of loops that 
detach during the evolution of the network in interconnections 
of open strings and smoothing of their small-scale structure. 
The network evolves perpetually as both long segments and 
open loops move and oscillate at relativistic velocities. 

Strings are believed to be a sub-dominant species in the 
matter-energy balance of the Universe. Analytical calcula- 
tions and numerical simulations indicate that string networks, 
soon in the course of their evolution, can reach a scaling be- 
havior where a typical distance between the strings increases 
in proportion to the horizon scale f3-TTo(] ■ This corre- 
sponds to the density of strings p decreasing with redshift z as 
p(z) oc d^ 2 (z) although one should keep in mind that these re- 
sults were obtained for radiation- and matter-dominated eras 
with no contribution from the vacuum energy. In the matter- 
dominated era this dependence coincides with that for the 
cold matter p(z) oc (1 + z) 3 , which appears to be a natural 
behavior for networks dominated by non-interacting loops, 
whose density would decrease solely due to the universal ex- 
pansion. For subsequent calculations, we will use both rela- 
tions, p(z) oc d^ 2 (z) and p(z) oc (1 + z) 3 ; the law followed by 
the strings in the actual universe is likely to be an interpola- 
tion between these two cases. 

Despite being a natural prediction of many cosmolo gica l 
theories, cosmic strings have not been observed yet 
which raises an obvious question about the origin of this 
discrepancy. Attempts to answer it would benefit from an 
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estimate of the actual density of stings in the real Uni- 
verse or, in the absence of their detection, an upper limit 
on this parameter. However, observational estimates of this 
kind are surprisingly scarce llii-jill]. In a recent paper [Tg| |. 
we constrained the local (at ~ 1 kpc scale) density of light 
(10 -16 < Gfi/c 2 < 10 -10 ) cosmic string loops using their 
observational signatures in pulsar timing and precision pho- 
tometric surveys. These constraints were made possible by 
significant enhancement in the local density of strings due to 
clustering of string loops expected in the Galaxy [2(j|. How- 
ever, this enhancement is subject to theoretical uncertainties 
that are hard to quantify at the current level of our under- 
standing of cosmic strings. 

In the present work, we derive robust observational upper 
limits on the average cosmological density of cosmic strings 
that are independent of any clustering effects and apply to 
both open segments and closed loops of strings. Instead of 
the local enhancement, we rely on giant distances to and a 
large number of extragalactic objects, namely quasars, used 
in deriving these constraints. Our method is based on the 
statistical analysis of quasar variability obtained for a large 
sample of quasars from the SDSS catalogue [2l[. Lensing by 
cosmic strings that are heavy enough (Gp,/c 2 > 10~ 14 — 10~ 12 , 
though this is somewhat quasar-model-dependent) would lead 
to an excess of twofold jumps in the distribution of brightness 
variation between two observational epochs. Hence, absence 
of any such features in the observed distribution allows us to 
infer robust upper limits on the density of strings. 

The paper is organized as follows. In Section [H] we elabo- 
rate on the idea above to see how the observed distribution of 
the variability in an ensemble of quasars can be used to put 
an upper limit on the probability of string lensing. Section Hill 
relates this limit to the density of strings by calculating the 
probability of lensing as a function of string and source pa- 
rameters. Then, in Section IIVI we use observational data to 
obtain actual limits on the density of strings from the SDSS 
data. Finally, in Section [V] we present our results in Figure [2] 
and conclude with a short discussion. 

All numerical calculations are made for a standard flat cos- 
mological model with a cosmological constant A = 1 — £2, cold 
matter density fl = 0.27 [2^] and the present-day Hubble con- 
stant Hq = 71 km ■ s^ 1 • Mpc -1 [23j |; we assume that strings 
do not contribute appreciably to the energy budget. 
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II. METHOD: OBSERVATIONAL LIMITS ON 
PROBABILITY OF FLUX DOUBLING 

Cosmic strings produce a distinctive pattern of lensing; for 
a point source inside a narrow strip along the string a second 
positive-parity image appears in a duplicate strip on the other 
side of the string [24l-[2a] . In the following, we assume that at 
any given moment the source is crossed by at most one string. 
This is a sensible approximation given that otherwise lensing 
by cosmic strings would be ubiquitous and would likely have 
been detected by now. A formal demonstration of the plau- 
sibility of this assumption relies on the fact that the strips' 
intersections cover a fractional area oc r 2 in projection, where 
t is the optical depth to string lensing. For cosmologically 
far-away sources (z s ~ 1), the string lensing optical depth, 
like that for point lenses, is of order the fraction of critical 
density fl s in the strings [ijj, which is unlikely to be greater 
than unity. 

When an extended source being crossed by the string can- 
not be resolved, the observer sees a flux increase given by the 
flux in the part of the source that is momentarily inside the 
strip. As the source moves into the strip, the flux gradually 
rises from a flat unlensed 'bottom' to some maximum value 
and then falls off to the same bottom as the original source 
disappears behind the strip while its duplicate image leaves 
the duplicate strip. The exact shape of the light curve de- 
pends on the brightness distribution in the source and the 
maximum is determined by the size of the source in rela- 
tion to the strip, which can only be guessed in the case of 
quasars. However, for small enough sources that fit into the 
strip completely, the maximum increase is exactly twofold, 
which is Amo = 2.5 lg 2 m 0.75 m in terms of stellar magni- 
tudes. Moreover, string lensing light curves of such sources 
possess a characteristic extended 'plateau' at this level, its 
width is given by the time it takes the source to traverse the 
strip. 

The distinctive shape of the light curve readily imprints 
itself in the distribution density of the magnification /i of a 
small source lensed by a cosmic string. This function consists 
of a certain smooth component at 1 < /i < 2 and a pair of 
5-functions that correspond to the unlensed case ^ = 1 and 
the maximally lensed case /i = 2; the latter events have a 
non-zero measure due to the extended nature of the corre- 
sponding 'bottom' and 'plateau' of the light curve. The dis- 
tribution density of the magnification allows one to calculate 
the density p(Am) of the magnitude jump Am in two obser- 
vations due to a change in the magnification factor between 
the corresponding epochs. This function is even because of 
the symmetry between the two epochs and consists of three 
5-functions at Am = and Am = ±Amo on top of a smooth 
component p(Am) for —Amo < Am < Amo: 

p(Am) = p(Am) + PS(Am qp Am ) + QS(Am). (1) 

For sources of unknown brightness profile, it is not possible 
to calculate the smooth component of this distribution; how- 
ever, the amplitudes of 5-functions can be calculated rather 
straightforwardly as shown in the next section. If lensing by 
cosmic strings is a rare phenomenon, which we will assume, 
2P is essentially the optical depth (see Eq. I14p to lensing by 
cosmic string, which is r -C 1. 

The ensemble variability studies approach the question of 
the variability of celestial objects by comparing the magni- 
tudes of a large number of individual sources observed at a 



few (two or more) epochs and presenting various statistical 
measures of the individual magnitude change in the ensem- 
ble - its mean, variance, distribution density, autocorrelation 
function and the like [2ll. [27l - [33 | . It is often 'ergodically' as- 
sumed that these measures reflect those of individual sources 
to an extent given by the size and homogeneity of the obser- 
vational sample and the time span of the variability survey. 
However, certain statistical measures of the ensemble vari- 
ability have a value on their own. Of these, the distribution 
density /(Am) of the observed magnitude change will be par- 
ticularly important for our study. 

If the sources crossed by cosmic strings were not variable, 
/(Am) would be a direct observational estimate of the under- 
lying density of lensing magnification p(Am). The quasars, on 
the contrary, are observed to vary at different magnitude and 
time scales (e.g., [3ll. [3^ - [3r3 | ) . Nevertheless, lensing by cosmic 
strings might still be apparent in /(Am) of these objects as 
±Amo inter-epoch changes will be overrepresented in the ob- 
served variability sample. This is the essence of our method 
to constrain the cosmological density of cosmic strings. 

Mathematically, this can be expressed as follows. The 
string-induced variability is extrinsic to quasars and there- 
fore the observed magnitude change distribution /(Am) is a 
convolution of those due to strings p(Am) and due to intrin- 
sic processes in quasars s(Am) ; for p(Am) given by (QJ the 
convolution equates to 

/(Am) = Qs(Am) + Ps(Am =p Amo)) + s(Am) (2) 

with s(Am) being the convolution of the intrinsic density 
s(Ara) and the (unknown) smooth component of lensing den- 
sity p( Am). Figure[T]shows an example of the observed distri- 
bution density /(Am) derived from the data presented in 
(see Section llVI for details); it shows little evidence for any 
excess at ±Amo, which will be used below to infer an upper 
limit on the density of cosmic strings. 

Equation ((2| immediately gives a handle on the parameter 
P related to the density of cosmic strings: 

p = /(Am) - Qs(Am) - s(Arn) 
s(Am + Am ) + s(Am - Am ) ' 

We do not know what the intrinsic variability s(Am) is and 
therefore cannot distill the string signal from the observed 
/(Am) directly. However, if we assume that all of the vari- 
ability at a certain level Am comes from strings, this clearly 
gives us an upper limit on their contribution to the variability, 
which can be used to infer robust constraints on the popula- 
tion of strings: 

P< fJ^l • (4 ) 

~ s(Ato + Amo) + s(Am - Amo) 

this inequality is valid irrespective to the assumptions on P 
because both neglected subtrahends in the numerator of the 
fraction in (J3} are non-negative. 

To deal with the denominator we assume that lensing by 
strings is rare; this is a sensible assumption as discussed 



There is also a distribution of observational uncertainties but 
it can be absorbed into the intrinsic variability distribution; we 
assume that for most quasars in the sample this distribution is 
only weakly dependent on their actual brightness. 
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FIG. 1: The observed distribution densities /(Am) of the 
brightness variation Am between photometric and spectro- 
scopic measurements by the Sloan Digital Sky Survey (SDSS) 
for > 25 000 quasars. The distributions in three SDSS pass- 
bands, i, r and g, are derived from the results of [2l|] using 
a procedure described in Section IIVI The error bars shown 
in the figure correspond to the 'Poissonian' square roots from 
the number of quasars in each Am bin. 



above. In this case the amplitude P <C 1, s(Am) <C s(Am), 
Q ~ 1, and the observed variability distribution density 
/(Am) is very close to the intrinsic one s(Am) - except, pos- 
sibly, at points Am = ±Amo, where a contribution Ps(0) 
due to an excess of Amo jumps from the lensing light curve 
plateau might be expected. It therefore makes sense to use (Q 
at one of those points to derive an upper limit P on the param- 
eter P. The denominator at these points can be approximated 
by the observed function /(Am) and one has 



/(±Am ) 



/(O) + /(±2Am ) 



/(±Am ) 
/(0) ' 



(5) 



the last step reflects the observational fact that / measured 
at Am = ±2Amo is orders of magnitude lower than at zero 
where it peaks (cf. Figure [1]). 

The constraints obtained in this way can be further re- 
fined and potentially even turned into assertive estimates for 
the string population properties by including additional in- 
formation such as dependence of the observed distribution of 
magnitude change on source parameters or inter-epoch time 
lag. This can be accomplished by calculating the probabili- 
ties of the observed data given model parameters and using 
the Bayes theorem to infer the reverse. However, such an en- 
deavor would inevitably require a model for the distribution of 
the intrinsic variability s(Am) and its dependance on source 
parameters, time lag or whatever else that is included in the 
analysis of the overall observed variability. In this study, we 
will use a simples approach outlined above, which is inde- 
pendent of authors' ignorance of the intrinsic variability of 
quasars though can only provide upper limits on the density 
of strings. 



III. MODEL: PROBABILITY AS A FUNCTION 
OF STRINGS POPULATION 

The amplitude P is the probability that the magnification 
Am jumps by Amo between the two observational epochs, t 
and t + At. Since Amo is the maximum brightness increase 
due to string lensing, the only configuration that corresponds 
to this jump is that where the source is completely inside 
the strip in one of the observations and completely outside 
the strip in the other. Because of the symmetry between the 
two epochs we can assume that it is the first observation when 
the source is inside the strip and the second when it is outside 
thereby replacing At with is absolute value: 



2P = V [Am(t) = Am and Am(i + |At|) = 0] 



(6) 



To estimate this value we first introduce the angular width A 
of the string lensing strip. According to [24|, Hg|, it depends 
on the tension fi of the string and its local inclination 6 to the 
line of sight : 



A = 8tt |sin<9 



G> As 
c 2 D os 



(7) 



where D oa and Di s are the (angular diameter) distances, 
respectively, from the observer and from the string to the 
source (along the line of sight); we use the average value of 
{ | sin 6*|) = 7r/4. It seems sufficient for our study to assume 
that the string segment responsible for lensing is long and 
straight compared to the angular size of the source; for a 
comprehensive study of lensing by general configurations of 
strings see [3?], HI] . 

Now let x be the initial epoch position of the projection of 
the source center onto the lens plane with respect to the strip 
median line (measured towards the outer edge of the strip such 
that the string itself is at x 3 = —A/2). The position of the 
source in the second epoch is then x + f3±c\At\/ (1 + z\)D i, 
where /3±c is the orthogonal (to the string) component of 
transverse (to the line of sight) velocity of the string w.r.t. the 
source; factor (l+zi) -1 corresponds to the dilation of observed 
time lag At from the lens plane at redshift Z\. Cosmic strings 
are exp ected to move relativistically, /3 ~ 0(1) [H; following 
[l9l.l39||. we use /3± = 0.3 in subsequent calculations 2 . 

Since we assume that the source is lensed by at most one 
string, conditions in the argument of (JfjJ) require that the cen- 
ter of the source is inside the strip by a margin of at least the 
source size r± — R±/D os (R± is its linear size) on the first 
observation and outside it by the same margin on the second 
epoch: 



< A/2 - r± 

h /3 ± c\At\/(l + z x )D ol > A/2 + r ± 



(8) 



Taken together, they restrict x to a narrow strip of width £, 
which is the lowest of A — 2r± and /3±c\ At\/ (1 + z\)D i — 2r± 
as long as this lowest is positive, and zero otherwise: 



£ = max < 0, min 



P±c\At\ 



' (1 + z\)D 



2r, 



(0) 



2 We assume that /3j_ > 0, i.e. the source is moving away from the 
string; this is not restrictive due to time symmetry mentioned. 
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The probability that a randomly placed source will lie within 
the strip of this width parallel to a string in an infinitesimally 
thin slice of string network with local number density pj[i is 



dr = ^DoidDoi = a 



3Hl 



£j(,zi)£-D idD 



(10) 



where D is the proper distance along the line of sight param- 
eterized by the slice redshift z\. 

In the formula above we also introduced the current cos- 
mological density of strings a and its dependence on redshift 
w(z\) such that the proper density p{z\) = cj(zi)f2 s 3//o/87rG. 



We use two models for oj(z\) 
solutions 



dh(0) 



and pressureless dust: 
oj(z) = (1 + zf. 



that corresponding to scaling 



(11) 



(12) 



As discussed in the Introduction, there is currently no consen- 
sus on the relative contribution of closed and open strings to 
the energy budget of the Universe, but whatever the contri- 
butions are, the final result for the density constraints can be 
obtained by interpolating those derived from the application 
of (ll]) and (fT2jl . 

With equations ((9l I10[) we can now write down the prob- 
ability 2P of a twofold magnification jump in any of the in- 
finitesimal slices. It is then given by the integrated optical 
depth r along the line of sight to the source 



T 

a 



8ttG> 



d-D(zi) oj(z{)D(zi)£ (21, z a , /i, R±, \At\) (13) 



according to 



2+^(r 2 ); 



(14) 



this probability is placed symmetrically in Am = ±Amo, 
hence the factor 1/2 in front of the brackets. 



IV. APPLICATION: OBSERVATIONAL LIMITS 
ON STRING DENSITY 

In order to put an upper limit Q a on the density of cosmic 
strings one now can simply equate the observational upper 
limit P on the probability of lensing to its model estimate P 
given by (| 14p in the limit r<l: 



a 



2P 

r/a 



(15) 



The numerator of the fraction above can be estimated us- 
ing from the distribution density /(Am) of quasar bright- 
ness variations, which can be calculated directly from the ob- 
servational data. In this regard, SDSS quasar survey [13] pro- 
vides an invaluable observational sample, where the bright- 
ness of tens of thousands of quasars is homogeneously mea- 
sured in a number of optical passbands and could be com- 
pared against an equally homogeneous sample of brightness 
estimates derived from the quasar spectra, which are obtained 
month and years after the photometric observations. 



TABLE I: Observational estimates for the distribution den- 
sity of quasar brightness variation at Am = 0, ±Amo in three 
SDSS passbands and corresponding estimates for P. For con- 
sistency, we calculate the central value of the distribution den- 
sity /(0 m ) = [/(-0.05 m ) + /(0.05 m )]/2; the value of P quoted 
in the table is the average between values corresponding to 
Am = — Amo and Am = Amo according to ([5]). 

Passband /(CP) /(-0.75 m ) /(0.75 m ) P 
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Such an analysis has indeed been done for N = 25710 SDSS 
quasars by Vanden Berk et al. [2l[ and we use their data on 
quasar variability in SDSS passbands g, r and i to derive con- 
straints on the string population. The authors of the cited 
study do not explicitly quote estimates on the probability 
distribution density /(Am) and we do not possess sufficient 
resource to re-reduce the publicly available SDSS data to de- 
rive /(Am) independently but it can be readily obtained from 
Figure 3 in the PDF version of [2l[. To do so, we manually 
counted the data points corresponding to individual quasar 
measurements in the scatter plots of the figure in 0.1 m -wide 
bins for |Amj > 0.4 m or summed the heights of the respective 
histograms for |Am| < 0.4 m 3 , and then divided them by the 
bin width and the total number of quasars in the ensemble. 

The distribution densities of quasar variability in three 
passbands obtained in this way are plotted in Figure [1] while 
Table U presents corresponding values /(Am) at points of our 
interest and an estimate for P in each passband. Since lens- 
ing is achromatic (as long as the string is heavy enough, such 
that A > 2r± for all passbands) and intrinsic variability is 
not, we are free to choose the lowest P to use in (|15f) : this is 
P = 3.2 ■ 10 -3 , which corresponds to passband r. 

The denominator of the fraction in (|15l) is given by the 
right-hand side of l)13|l . which depends on the source red- 
shift, and we therefore need to take an average with respect 
to the distribution of the observed quasars. The quasar sam- 
ple of [2jl includes most of the quasars in SDSS Data Release 
1 Quasar Catalogue (4]| and a substantial fraction of quasars 
observed by SDSS that were not included in SDSS DR1. The 
properties of individual quasars in the entire sample used in 
that study do not appear to have ever been detailed in a pub- 
lication and therefore we use SDSS DR3 QSO catalogue [i^l 
as a proxy to the statistical properties of the true sample. 
We have verified numerically, that our results do not change 
significantly if we use either DR1 or DR5 catalogues in 
averaging r/a (numbers do get higher when using more re- 



3 The quality of the figure does not allow us to use a consistent 
counting approach in the entire domain of Am - inner regions 
(|Am < 0.4 m ) of scatter plots suffer from considerable confusion 
of data points while the linear scale of the histograms makes 
them hardly readable for |Am| > 0.6" 1 . However, where this 
comparison is possible, at bins centered at ±0.45 m and ±0.55 m , 
the numbers agree to within a few percent, which is acceptable 
given somewhat low-tech approach employed in the absence of 
published digital data. 
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cent, deeper versions of the catalogue but only by 2-3 per 
cent). 

Another source of uncertainty is the physical size 2R± that 
produces most of the observed flux. The size of the quasar 
affects our results significantly by limiting the sensitivity of 
our estimate as a function of string tension /i. Estimates 
on these quantities vary appreciably in the literature, de- 
pending on the method used; reverberation mapping seems 
to favor sizes in the range R ~ (fO 16 - 10 17 )cm 
while microlensing tech niq ues give somewhat smaller values 
of R ~ (f0 15 -f0 16 ) cm [49M53t| - Both methods are not model- 
independent and the estimates are expected to correlate with 
individual properties of QSOs such as the luminosity or the 
mass of the central black hole. In the apparent absence 
of a better option, we perform our calculations using three 
representative values of 2R ± € {I0 15 , 10 16 , 10 17 } cm treating 
10 16 cm as a fiducial estimate. 

Finally, the value of the time lag At between two obser- 
vational epochs in the observer frame cannot be read from 
the results of [U directly, which also introduces some uncer- 
tainty. The value of At affects our results directly via Q and 
therefore need to be fixed to perform calculations. From the 
visual inspection of Figure 4 of [2lj] it is clear that typical time 
lag in the source frame At/(1 + z s ) ~ (100 - 200) days. The 
average redshift of SDSS quasars is z s » 1.5. We therefore 
take At = 150 days x 2.5 w 3.2 ■ 10 7 s. 



V. RESULTS AND DISCUSSION 




FIG. 2: Upper limit on the average present-day cosmologi- 
cal density of light cosmic strings as a function of its tension 
fi based on the quasar variability distribution shown in Fig- 
ure [T] Three pairs of curves are shown for the assumed source 
size 2R± of 10 16 cm (black, middle), 10 17 cm (blue, top) and 
10 cm (red, bottom). Solid lines correspond to 'scaling' evo- 
lution of string density with redshift according to {TTJ, dashed 
ones assume a 'dust-like' law (|12[) ; they do not differ much. 
A thin grey line shows the constraints obtained in [li| from 
the local effects of strings. 



Figure [2] presents the upper limits on the average cosmo- 
logical density of strings set by the statistics of the observed 
variability of more than 25 thousand SDSS quasars as a func- 
tion of string tension Gfi/c 2 . Depending on the assumed size 
of the source, which turns to be the major parameter of the 
method developed in this paper, the region where the strings 
density is usefully constrained extends for up to five orders of 
magnitude. At the same time, our constraints are only weakly 
dependent on the assumed behavior of string density with red- 
shift, mostly because prescriptions (|11[) and (|12p started to 
diverge from each other relatively recently, when the vacuum 
density began to dominate in the Universe. 

The most stringent constraints are obtained at G/i/c 2 ~ 
(10 _13 -10 -11 ) where the upper limits reach below the level of 
Sl s = 0.01. We note that these are rather weak limits for open 
topological cosmic strings because their cosmological density 
is believed to be of order Q, B ~ (10 — 100)G/x/c 2 based on 
the reconnection argument leading to scaling solutions (e.g., 
0). However, this argument might not be applicable to fun- 
damental strings, for which the reconnection probability is 
highly model dependent and can be significantly lower than 
unity [54ll5"5l|; neither it is clear how it applies to cosmic string 
loops, which are essentially a product of the mechanism that 
ensures that the density of open string has a scaling behavior. 
Our observational upper limits apply to all kinds of strings, 
whether topological or fundamental, open or otherwise, and 
therefore provide strong constraints in the case of string loops 
and fundamental strings, independent of theoretical uncer- 
tainty. Moreover, our results might potentially be useful for 
constraining the unknown reconnection probability of funda- 
mental strings. 

The upper limits obtained here rely on an extensive dataset 
of SDSS quasars and are robust in the sense that they do not 



depend on any local enhancement of the density expected for 
string loops in the Galaxy [2(|] or a particular model of quasar 
emission, as long as the source remains sufficiently small com- 
pared to the lensing strip width. This method is also rather 
insensitive to the photometric accuracy of observations be- 
cause 0.75 m -wide jumps in brightness are fairly obvious by 
any standard. And if follow-up observations could be per- 
formed when a survey telescope sees a sudden twofold increase 
in the brightness of a quasar, it might even be possible to con- 
firm or rule out string nature of this increase to a high level 
of confidence. 

However, all this robustness also means that upper limits 
obtained in this paper could not be improved much by simply 
increasing the number of data points in the sample, which 
is expected to grow by orders of magnitude with the launch 
of next generation photometric surveys, such as LSST [56( [. 
Rather, it is the study of quasar intrinsic variability that can 
bring most benefit for this cause. The variation distribution 
density presented in Figure [1] does not appear to show any 
notable features around Am = ±0.75 m which would make 
these points stand our from the smooth decrease of /(Am) 
with |Am|. Therefore, the upper limits derived in this paper 
are most likely a significant over-estimation of the true den- 
sity of strings because what is conservatively interpreted as a 
string 'signal' here is most likely the quasar variability 'noise'. 
It will not simply walk away if one takes more quasars. To 
make further progress here, one needs to learn to separate 
two contributions - e.g., using the varying dependence of the 
two effects on various parameters, such as the source redshift, 
color, luminosity and so on. We can easily calculate how 
the statistics of string lensing effect should depend on these 
parameters but at present are less certain when it comes to 
quasar intrinsic variability. 
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Nevertheless, the science of quasar variability is advancing 
fast and our understanding of it might soon be sufficient for 
distilling the string effect from the observational data. This 
approach can also be followed in the analysis of future datasets 
from the GAIA mission [57| on the variability of stars, which 
are less distant but photometrically more stable and much 
more numerous than SDSS quasars. Moreover, at present we 
do not know the density of cosmic strings in the Universe and 
cannot predict just when the twin peaks of strings signal at 
Am = ±Amo will begin to show up in the growing datasets 
of ensemble variability studies. 
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